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Abstract — The paper presents a prediction-identification model 
based adaptive control method for uncertain systems with time 
varying parameters in the presence of bounded external distur- 
bances. The method guarantees desired tracking performance for 
the system’s state and input signals. This is achieved by feeding 
back the state prediction error to the identification model. It 
is shown that the desired closed-loop properties are obtained 
with fast adaptation when the error feedback gain is selected 
proportional to the square root of the adaptation rate. The 
theoretical findings are confirmed via a simulation example. 

Index Terms — Time-varying parameters, guaranteed transient 
performance, disturbance rejection; 

I. Introduction 

Adaptive control has been considered as a promising tech- 
nology to improve stability and performance of uncertain 
systems. However, limitations of conventional adaptive meth- 
ods (see for example [1]) have prevented them to be widely 
adopted in safety-critical systems. 

During past two decades majority of the efforts have been 
directed to improving the transient of the tracking error 
(see for example recent results in [2], [4], [5], [9] for the 
systems with time varying uncertainties), but not the control 
signal, the behavior of which significantly contributes to the 
aforementioned limitations. 

These limitations have been addressed in the C\ adaptive 
control framework [3], It has been shown that the desired 
transient can be obtained via fast adaptation and a low-pass 
filter, which a priori sets the bandwidth, within which the 
uncertainties in the system can be compensated for. 

An alternative method, which guarantees desired transient 
behavior of the closed-loop system, has been proposed in [8]. 
It is based on the modification of the reference model by 
the tracking error feedback, and is called modified reference 
model MRAC (M-MRAC). The idea behind the method was 
to drive the reference model toward the system proportional 
to the tracking error, thus preventing the system’s attempt 
to aggressively maneuver toward the reference model in the 
initial stage of the process. It turns out that the error feedback 
gain determines the damping in the control signal dynamics, 
whereas the adaptation rate determines the frequency. This 
allows the designer to choose proper values for the parameters. 
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In this paper, we extend the approach to the systems with 
time varying parameters using a prediction (or identification) 
model based approach. The parameters are assumed to be 
bounded and vary with bounded derivatives, no matter how 
large these bounds are. The parameter estimates are generated 
using the state prediction error as in the case of conventional 
indirect adaptive control schemes, which is the reason to name 
the approach indirect M-MRAC. However, our prediction 
model differs from the conventional ones by a prediction 
error feedback term, which turns out to play the same role 
as the tracking error feedback term plays in the direct M- 
MRAC approach. Hence, the desired closed-loop behavior 
can be achieved with fast adaptation by selecting a proper 
error feedback gain, which also separates the time scale of the 
adaptive estimation from that of the system’s dynamics. 

The rest of the paper outlines the properties of the proposed 
indirect M-MRAC control architecture and demonstrates the 
benefits of it in a simulation example. 

II. Problem Statement 
Consider the system 

x{t) = Ax(t) + BA(t) + K(t)g(x(t)) + d(t)} (1) 

with a:(0) = xq, where x £ R n is the state of the system, 
u £ R q is the control, A £ R nxn and B £ R nxq are 
given constant matrices with A being Hurwitz and (A,B) 
controllable, g : R" -A R p is continuously differentiable, 
A : R+ -£ R qxq is positive definite with bounded and 
piecewise continuous unknown entries, which have bounded 
derivatives, K : R + -A R qxp is an unknown parameter matrix 
with bounded and piecewise continuous entries, which have 
bounded derivatives, and d : R 1 — > R q is bounded and 
piecewise continuous disturbance with a bounded derivative. 

The control objective is to design a control input such that 
the system (1) tracks the reference model. 

x m (t) = Ax m (t) + Br{t) (2) 

with x m (0) = x m o, where r : R + -A R q is bounded and 
piecewise continuous command with a bounded derivative. 
We notice that the system (1) can be represented in the form 

x(t) = Ax(t) + Br(t) 

+ BA(t)[u(t) + @(t)f(x,r) + d(t)] , (3) 

where 0(f)/(x,r) = K(t)g(x(t)) — A ~ 1 (t)r(t). Let 

l|A(f)|ko<A*, ||0(f)|Uoo < 0*. H*)ll c„<d* 

l|A(t)|| £oo < as, ||6(t)| | £oo < r d , Hi)|| £oo < d* d (4) 



III. Prediction Model 


We introduce the following adaptive prediction model 

x(t) = Ax(t) + Br(t) (5) 

+ BA(t)[u(t) + Q(t)f(x, r) + d(t)\ + kx(t) 


with £c( 0 ) = Xo, where x(t) = x(t) — x(t) is the prediction 
error, k > 0 is a design parameter, A(f), 0 (f) and d(t) are 
the estimates of the unknown quantities, generated according 
to adaptive laws 

0(f) = 7 Pr (©(£), B T Px(t)f T (x, r )^ 

A(f) = 7 Pr (k(t), B T Px(t)[u(t) + 0(f)/(x,r)] T ) 

d(t) = 7 Pr^d(f), B T Px(t)^j , ( 6 ) 

where 7 > 0 is the adaptation rate, P = P T > 0 is the 
solution of the Lyapunov equation A T P+PA = — Q for some 
Q = Q t > 0, and Pr (-, •) denotes the projection operator [7], 
which is defined as Pr (9,y) = [I — G(9)\y, where 

! 0 , if <p{0) < 0 

0, if^)>0, V^ T (0)y<O 

if ^ ] - °’ > 0 

with the notation V<p(0) = , and the smooth convex 

functions ip(9) is given by ip(9) = gl 9max with 0 m ax 
denoting the norm bound imposed on the parameter matrix 
9 and eg denoting the convergence tolerance. The projection 
operator has the following properties 
Lemma 3.1: [7] Let 9 0 € = {9 G R" | ip{6) < 0}, and 

let the parameter 0 (f) evolve according to the dynamics 

9(t) = Pr(9(t), y), 0(£ o )efi. (7) 

Then 1) 0(f) e ffi = {9 £ R” | ip(9) < 1} or ||0(£)|| < 
9* for all t > to, where 9* = /I + eg 0 max , 2) [0(f) — 
0 o ] T [Pr(0(f), y)-y]< 0 for all f > f 0 . 

It is straightforward to verify that 


x(t) = (A — kl)x(t) + BA(t)[Q(t)f(x, r) + d(f)] 

+ BA(t)[u(t) + Q(t)f(x,r) + d(t)], ( 8 ) 

where 0(f) = 0(f) — 0(f), A(f) = A(f) — A(f) and d(f) = 
d(f) — d(f) are the estimation errors. 

Lemma 3.2: If Xq = xq, then the prediction error x(t) 
satisfies the bound 

W' (9) 

where c = c i + f|, c\ = 4A *d* 2 + 4A *$* 2 + 4A* 2 , and 

c 2 = 4A*i?*<?5 + 4X*d*d d + 4X d d* 2 + ±X* d d* 2 . 

Proof: The derivative of the candidate Lyapunov function 

V(t) = x T (t)Px(t) + 7” 1 <f T (f)A(f)d(f) ( 10 ) 

+ 7 _1 tr ^0 T (£)A(£)0(f) + A T (f)A(f)j , 


computed along the trajectories of the prediction error dynam- 
ics ( 8 ) and the adaptive laws ( 6 ), satisfies the inequality 

V(t) < —x T (t)Qx(t) — 2kx T (t)Px(t) (11) 

+ 2 7 _ 1 tr(0 T (£)A(£)0(£)) + 2 7 _ 1 d T (f)A(f)d(f) 
+ 7 ^ 1 d T (f)A(f)d(f) + 7 _1 tr ^0 T (f)A(f)0(f)^ . 

The projection operator in the adaptive laws ( 6 ) guarantees 
the inequalities || 0 (f)|| < d*, ||A(f)|| < A* and ||d(£)|| < 
d*. Therefore || 0 (£)|| < 20*. ||A(t)|| < 2 A*,Jd(f)|| < 2d* 
and 2 tr( 0 T (f)A(f) 0 (f))+ 2 d (f)A(f)d(£)+d (f)A(f)d(f) + 
tr ^ 0 T (f)A(f) 0 (f)^ < c 2 . That is 

V(t ) < — x T (t)Qx(t) — 2kx T (t)Px(t) + 7 _ 1 c 2 . (12) 

On the other hand we have d (f)A(f)d(f) + 

tr ^0 T (f)A(f)0(f) + A T (f)A(f)^ < ci. It follows that 

V(t) < x T (t)Px T (t) + 7 _ 1 ci. Therefore, if V(t) > 7 _ 1 c, 
for some r then x T (t)Px t (t) > ( 2k r y)~ 1 C2 , which 
implies that V(t) < 0. Since a;(0) = 0 it follows that 
V(0) < 7 _1 ci < 7 _ 1 c. Therefore V(t) < 7 ~ 1 c for all £ > 0. 

Since ||*(f )|| 2 < x T (t)Px{t)/\ min {P) < V(t)/X min (P), 
the inequality (9) follows. ■ 

It can be observed from Lemma 3.2 that the state prediction 
error can be decreased as desired by increasing the adaptation 
rate 7, when the prediction model is precisely initialized. The 
next lemma shows that the initialization error results in an 
additive exponentially decaying term. 

Lemma 3.3: If Xq ^ then x(t) satisfies the bound 

+ (13) 

where C 3 = |V(0) — ^|, and V(t) is defined by (10). 

Proof: Using the same V(t) as in Lemma 3.2 and 
following the same steps one can arrive to the inequality 

V (f) < -2 k[V (f) - 7 _ 1 ci] + 7 _1 c 2 , (14) 


integration of which results in 

c 


V(t)< 


v ^) - - 


— 2 kt 1 C / —2 kt 1 t' /in 

e -I — < C 3 e -I — . (15) 

7 7 


Recalling that ||S(f)|| 2 < U(f)/A m i n (P), we readily obtain 
||*(i)|| < \l~ 7777 \l c 3 e 


Amin (P) 


—kt 


7 


(16) 


Taking into account the inequality \J a + b < ,/ci + \/b for any 
a > 0, b > 0, the bound (13) is concluded. ■ 

Since the effect of the prediction model initialization error 
decays exponentially with the rate k, which is assumed to 
be set to large values for the fast adaptation, in the next 
derivations we assume that xq = x () . 



IV. Control Design 


Sideslip angle in degrees 


Since the reference model is designed to satisfy the ro- 
bustness and performance specifications, one would naturally 
select the control signal 

u(t) = - d(t) . (17) 

to achieve the control objective, if the system (3) were com- 
pletely known. Obviously, u(t) is not implementable, therefore 
its adaptive version 

u{t) = —Q T (t)f(x, r ) - d(t) (18) 

is used. When this control signal is applied, the prediction 
model (5) reduces to the modified reference model introduced 
in the M-MRAC architecture, that is 

x(t) = Ax(t) + Br(t) + kx(t) , (19) 



Bank angle in degrees 



Fig. 1. M-MRAC angle tracking performance with 7 = 1000. 


Lemma 4.1: Let the system (3) be controlled by the adap- 
tive control u(t) (18). Then closed loop signals are bounded. 

Proof: Under the action of the adaptive control u(t), 
the error between the prediction model and reference model 
emit) = x(t) — x m (t) satisfies the equation 

e m (t) = Ae m {t) + Xx{t ) , (20) 

Since A is Hurwitz, and x(t) is bounded according to Lemma 
3.2, it follows that e m (f) is bounded. Since the input r(i) 
is bounded, the reference model’s state x m (t) is bounded, 
therefore the predicted state x(t) is bounded as well. Then, it 
follows that the system’s state x(t) is bounded. The parameter 
estimates are guaranteed to be bounded by the projection 
operator, therefore tt(f) is also bounded. ■ 

Lemma 4.2: Let the system (3) be controlled by the con- 
troller (18), which is defined by the prediction model (5) and 
the adaptive law (6). Then 

11*0)11 < /3ie -yit +/3 2 7"C (21) 

where u(t) = u(t) — u(t), and fa and fa are positive 
constants to be specified in the proof. 

Proof: It is easy to show that u(t) satisfies the equation 


’ uit) ' 



0 

I qxq 

' u(t) ' 

uit) 


. -7 Fit) Lit) 

qxq 



+ 

' 0 ' 
. 7 . 

21(f) + 

' 1 ' 
. k 

Z 2 it), 


where we denote p(t) = f T (x, r)f{x, r) + 1, F(t) = 
p(t)I qxq -H(t), H(t) = G(0)f T (x,r)f(x,r) + G(d ), L = 
B T PBA(t) (Lit) is positive definite), zi(f) = [pf)BQ P + 
p{t)Bj PAm\x{t), and z- 2 it) = — 0 T (f)/(:c, r) — d{t). Since 
x(t) is bounded, it follows that pit) is bounded. That is, there 
exists a positive constant a\ such that 1 < pit) < a\ for 
all 0 < t < oo. On the other hand, it follows from the 
definition of the projection operator that ||G(0)|| < 1 and 
||G(d)|| < 1. Therefore Fit) is bounded. Further, it follows 
from the dynamics (3) that x(t) is bounded. Therefore pit) 
and z%{t) are bounded. That is, there exist positive constants 


Yaw rate in degrees/second 





Fig. 2. M-MRAC rate tracking performance with 7 = 1000. 


a .2 and <23 such that \pit)\ < 02 and ||z 2 (i)|| < 0:3 for all 
0 < t < 00 . Hence, (22) can be considered as a second order 
linear system in u(t) with time varying coefficients, where 
the adaptation rate 7 determines the frequency of u(t) and the 
feedback gain k determines the damping ratio. We notice that 
selection of the initial parameter estimates inside the convex 
sets defined by the projection operator results in fT(f) = 0 
on some initial interval [0 fi]. Therefore, F(t) = pit) \ q on 
[0 t\]. Let ao = aiA 9 +V , where A 0 = max t >o X(L(t)) and 
Ao = min t > 0 A(L(f)). Denoting E{t) = aoI g — p{t)L{t), we 
can write 


+ 


uf) 

u(t) 

' 0 ' 

7 


Zl(t) + 


0 

7 a oIqX(j 
1 
k 


u-qxq 

-kl 


qxq 


u(t) 

u(t) 


(23) 


22 it) +7 


WflXfl 

m 


uit) , 


the solution of which has an equivalent integral form 


uit) = ip it) [*(0) *(0)] T +7/0 ^ 2 (i- t)[z!(t) + 
Eir)uiT)\dT + J*[f>iit - t) + kfaf - r)]z 2 (r) dr, (24) 


Aileron deflection in degrees 


where = [ipi (t) ip 2 (t)] is the first row of the transition 
matrix of the LTI part of system (23). Following [8], we select 

k = 2VTOo ; (25) 

which results in the minimum norm ||'i/ , 2(f : )||£i = (700) _1 . 
For the same k, we have ||'0i(f) + At^2(i)||.Ci < 

4(7ao) _1//2 . Since ||.E(i)ll£oo = ao — A 0 , we obtain 1 — 
llllMOIkJ-KWIUoo = 2A 0 a^ 1 . Then, according to [8], it 
follows from the expression (24) that 

ll*(*)ll < 6i(||*(0)|| + ||ix(0)||) e -^ 

+ (26) 

where b\ is a positive constant and 

v\ = — 2^ (v^- \J a$ — Ao^j . (27) 

From the definition of 27(f) and Lemma 3.2 we have 

IktWIUo. < (a 2 \\B T P\\+a 1 \\B T PA\\) ] l^^. 

Then, it is straightforward to obtain the bound (21) with 

„ Ml«0)|| + ||6(0)||) 

A = 2Ao 

- 4 «3\/ao + (q 2 ||^ t -P|| + ai\\B T PA\\)^c 
2Ao 2AoVAmi „(P) 

This concludes the proof. ■ 


V. Tracking Error 

In this section we derive a norm bound for the tracking error 
e(t) = x(t) — x m (t), which is given by the following lemma. 

Lemma 5.1: Let the system (1) be controlled by the con- 
troller (18), which is defined by the prediction model (5) and 
the adaptive law (6). Then 

||e(f)|| < ke~ ut + /3 4 7“5 , (28) 


/? 3 and /?4 are positive constants to be specified in the proof. 

Proof: It is straightforward to obtain the tracking error 
dynamics in the form 

e(f) = Ae{t) — BA(t)u(t) . (29) 


Since A is Hurwitz, it follows that there exist positive con- 
stants 62 and such that ||e At || < Therefore the 

following bound can be obtained 


||e(t)|| < fo 2 ||e(0)||e— ^ + 6 2 1| ^A(i)|| ^ f* e~^-^ ■ 
[/) ie~ VlT + /3 2 ^f~^]dT < & 2 ||e(0)||e _! ' 2t + b 2 \\BA(t)\\c oa 


_ e~ vit ) + % (1 - e-^h-i 


(30) 


which can be expressed in the form of (28) with 

(h = ^ 2 1| e(0) || + — r ||SA(f)||£ 

\v 2 - vi\ 

Pa = -b 2 fa\\BA{t)\\ Cao . 

Vl 
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Fig. 3. M-MRAC control signal time history with 7 = 1000. 
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Fig. 4. MRAC angle tracking performance with 7 = 1000. 


and v = min[i/i, v 2 \. ■ 

Remark 5.1: The proposed method guarantees the regula- 
tion of all error signals by increasing the adaptation rate, which 
is only subject to available computational power. Therefore, 
with fast adaptation the control objective is achieved without 
generating unwanted excursions and oscillations in adaptive 
signals. The effects of the external disturbances and parameter 
variations are compensated for by the fast adaptation, and the 
effects of the initial conditions decay exponentially. □ 

Remark 5.2: It can be observed that the dynamics of the 
reference model, the operating system and the tracking error 
have the same time scales determined by the matrix A. Hence, 
the reference model initialization error generates an additive 
exponential term & 2 ||e(0)||e _l/t in the bound of the tracking 
error with a rate of decay defined by the time constant of 
the reference model, since the adaptation process is much 
faster. The time scale of the prediction error dynamics (8) 
is determined by k, which is proportional to 77, and the time 
scale of the adaptive estimates is determined by 7. Therefore, 
for large values of 7 the time scale of the adaptive estimation 
process is separated from the time scale of the underlying 


closed-loop dynamics, which is not achievable by conventional 
adaptive methods [1], □ 

VI. Illustrative Example 

In this section, the advantages of the proposed indirect 
M-MRAC architecture are demonstrated in simulations for a 
dynamic model that represents the lateral-directional motion 
of a generic transport aircraft (GTM) [6], The nominal model 
is the linearized lateral-directional dynamics of GTM at the 
altitude of 30, 000 ft and speed of 0.8 M and is given by 

x(t) = A n x(t ) + B n u(t) , (31) 

where x = [/3 r p 0] T is the lateral-directional state vector, 
in which (3 is the sideslip angle, r is the yaw rate, p is the 
roll rate, <f> is the bank angle, and u = [<5„ <5 r ] 1 is the control 
signal that includes the aileron deflection S a and the rudder 
deflection 8 r , and the numerical values for A „ and B n are 
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Fig. 5. MRAC rate tracking performance with 7 = 1000. 
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The reference model is selected from the perspective of 
improving the performance characteristics of the nominal dy- 
namics and is given by the equation (2), where A = A n —B n K 
and B = B n N, with the feedback and feedforward matrices 
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Fig. 6 . MRAC control signal time history with 7 = 1000. 


The reference model is driven by a command, which is 
chosen to be a series of coordinated turn maneuvers. That 
is, sideslip angle command is set to zero and the bank angle 
command is chosen to be a square wave of the amplitude of 
15 degrees and of the frequency -f- rad/sec, which is filtered 
through a first order stable filter ^ |'| 0 . 

The uncertain model of GTM roughly corresponds to 28% 
loss of left wing tip at t = 0 sec, and 55% loss of rudder 
surface and vertical tail at t = 20 sec. Its dynamics are in 
the form of the equation (3) with piecewise constant 0(f) and 
A(f), and fix) = x. The corresponding numerical values are 


0(f) = 


-0.1820 0.0149 -0.1049 0 

0.0807 -0.0109 0.0168 0 

-0.2268 0.0209 -0.1053 0 
-0.8514 0.0692 0.0003 0 


f < 20 


f > 20 


A(f) = 



0.5401 

0.0167 ' 


-0.0632 

1.0524 


' 0.5413 

-0.0492 ' 


0.0408 

0.4225 


f < 20 


f > 20 


The external disturbance is chosen to be a sinusoid of am- 
plitude 0.1 and frequency 2n/3 rad/sec in the yaw channel and 
a square wave of amplitude 0.15 and frequency 7r/3 rad/sec 
in the roll channel. The disturbance magnitude corresponds to 
8.6 degrees of aileron deflection and 5.7 degrees of rudder 
deflection. In the definition of the projection operator the 
conservative bounds A* = d* = d* = 10 are used. 

First, a simulation is run with 7 = 1000, Q = I4 and k 
is computed according to (25), where we used conservative 
bounds Ao = 0.2 and A 0 = 2. Figures 1 and 2 display 
the tracking performance of the states. Clearly good tracking 
is achieved with the chosen controller gains, for which the 
control time history is presented in Figure 3. It can be observed 
that the adaptive control signal closely follows the reference 
signal given by (17). Small spikes in the control signal are 
attributed to the discontinuities of the disturbance. For the 
comparison purposes we also present the conventional MRAC 
performance with the same setup. It can be observed from the 
Figure 4 that MRAC achieves output tracking with small os- 
cillations in sideslip angle. However, the rates and the control 
surface deflection commands are experiencing unacceptable 
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Fig. 7. M-MRAC output tracking with 7 = 10000. 
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Fig. 8. M-MRAC input tracking with 7 = 10000. 

oscillations (see Figures 5 and 6 ). 

Next we increase the adaptation rate 10 fold. As it can be 
viewed form Figures 7 and 8 , the output and input tracking 
performances are substantially improved as predicted. Com- 
putations show that the tracking error is decreased more than 
v / T(l fold, implying the the derived bounds are conservative. 
Farther increase of adaptation rate to 7 = 100000 further 
improves the system’s input and output performance (see 
Figures 9 and 10), which verifies the theoretical derivations. 

VII. Conclusions 

We have presented indirect modified reference model 
MRAC (M-MRAC) approach to uncertain systems with time 
varying parameters and bounded external disturbances without 
imposing ’’slow variation” restriction on the system’s param- 
eters. The method uses a prediction error feedback term to 
speed up the adaptive estimation process, which results in 
predictable transient behavior for both state and input variables 
of the system. It has been shown that the unwanted high 
frequency effects of the fast adaptation in the control signal 
can be regulated by the proper choice of the feedback gain. 
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Fig. 9. M-MRAC output tracking with 7 = 100000. 
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Fig. 10. M-MRAC input tracking 7 = 100000. 
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